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摘 要 : 为 研究 气候 变化 背景 下 天 水 麦 积 山 风景 名 胜 区 1980 一 2019 年 油 松 径 向 生长 与 气候 因子 响应 模式 的 变化 
及 机 制 ,利用 树 轮 气候 学 方法 和 Vaganov-Shashkin 模 型 ,研究 了 气温 突变 前 后 天 水 麦 积 山 油 松 的 气候 响应 和 生长 过 
程 。 结 果 表 明 :(1) 麦 积 山 气温 于 1997 年 发 生 突变 ,突变 后 气温 显著 高 于 突变 前 。1980 一 1997 年 和 1998 一 2019 年 
油 松 与 气候 因子 响应 变化 呈 相 关 性 下 降 (5 月 气温 )、 上 升 (10 月 气温 .降水 ,7 月 降水 ,12 月 气温 ) 和 震荡 (6 月 降水 7 
月 气温 )3 种 模式 .(2) 气候 显著 变 暖 使 春秋 季 油 松 生 长 热量 供应 更 加 充足 ,生长 季 显 车 延长 ;夏季 油 松 生长 受到 高 
温 影 响 和 水 分 胁迫 。(3 ) 油 松 生长 与 气候 因子 响应 模式 的 变化 主要 是 气候 变 暖和 生长 季 的 变化 引起 ,如 果 气 温 持续 
上 升 其 响应 模式 可 能 进一步 变化 ,类 似 的 变化 过 程 可 能 也 在 其 他 油 松 生长 区 存在 。 
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北半球 中 高 纬度 地 区 树木 生长 受到 气候 因子 
的 强烈 影响 ,其 年 轮 中 记录 了 大 量 气候 信息 "”。 这 
些 信息 定年 准确 分辨 率 高 . 复 本 易 得 ,使 树 轮 在 树 
木 生长 对 气候 响应 和 区 域 气候 重建 等 研究 中 得 到 
EAH”, 树木 年 轮 学 研究 中 往往 认为 区 域 影响 
树木 生长 的 因子 是 稳定 的 ? ,通过 研究 主要 限制 
子 与 树木 生长 的 关系 可 以 建立 两 者 间 线 性 转换 函 
BCU ,并 基于 此 重建 气候 。 但 是 ,20 世 纪 中 叶 以 来 ， 
北半球 高 纬度 地 区 许多 树木 年 轮 记 录 中 发 现 了 被 
称 为 “分 离 效 应 ”的 异常 现象 , 即 树木 生长 对 气温 的 
敏感 性 发 生 了 异常 下 降 *” ,这 一 现象 说 明 树 木 生长 
与 气候 关系 不 一 定 是 稳定 的 。 青 藏 高 原 ”\ 川 西 臣 
龙山 ”秦岭 牛 背 深 “的 研究 也 发 现 了 ”分离 效应 ” 
的 存在 ,证明 这 一 现象 在 中 低 纬 度 和 高 海拔 区 域 也 
有 体现 。Babst 等 ”基于 全 球 2710 个 树 轮 样 点 的 研 
究 更 是 指出 1930 一 1960 年 和 1960 一 1990 年 干 冷 地 
区 树木 生长 的 气温 响应 显著 降低 ,而 几乎 全 球 树木 
生长 都 受到 了 更 强烈 的 水 分 胁迫 ,持续 的 气候 变化 
可 能 导致 全 球 树木 生长 的 气候 驱动 因素 的 再 分 
配 。 因 此 ,深入 理解 气候 变化 背景 下 树木 生长 与 气 
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候 响 应 关系 必须 考虑 多 因子 , 非 线 性 的 作用 过 程 ””， 
采用 树 轮 生理 模型 是 一 个 可 行 的 解决 方法 "1。 

Vaganov-Shashkin 模 型 (VS 模型 ) 基 于 逐日 气 
温 、 降 水 和 太阳 辐 照 度 模拟 季节 性 形成 层 活动 和 细 
胞 生成 过 程 " ,得 出 树木 逐日 生长 速率 .生长 主要 
限制 因子 和 生长 季 日 期 等 树木 生长 的 重要 参数 ,并 
可 以 对 比 实际 采集 的 树 轮 样本 验证 其 可 靠 性 ， 
此 ,在 探讨 树木 生长 过 程 和 影响 因素 的 研究 中 得 到 
了 广泛 应 用 。 目 前 ,国内 外 基于 VS 模型 已 经 开展 了 
大 量 的 研究 ,例如 ,Chen 等 应 用 VS 模型 指出 气温 
导致 的 干旱 是 加 拿 大 西部 白云 杉 生长 的 主要 限制 
因子 , 陈 兰 等 "运用 VS 模型 探究 了 秦岭 太白 红 杉 径 
向 生长 对 气候 因子 的 响应 ,发 现 生 长 季 (4 一 8 月 ) 的 
气温 .生长 初期 及 7 月 .8 月 的 降水 是 限制 秦岭 太白 
红 杉 生长 的 主要 气候 因子 。VS 模 型 可 以 模拟 树木 
生长 的 逐日 特征 和 各 因素 影响 ,运用 VS 模型 有 望 解 
释 树木 生长 与 气候 因子 响应 关系 的 变化 机 制 , 但 目 
前 相关 人 研究 还 较 少 。 

本 文选 择 位 于 西 秦 岭 北 坡 , 兼 受 东南 季风 E 
南 季风 和 高 原 季 风 影 响 ,气候 系统 独特 复杂 、 对 气 


E 与 资源 环境 . E-mail: yaodaijun@stumail.nwu.edu.cn 


http: //azr.xjegi.com 


202303.00121v1 


chinaXiv 


20 第 大 


候 变 化 异常 敏感 的 天 水 麦 积 山地 区 作为 研究 区 ， 
以 常用 于 树木 年 轮 学 研究 的 油 松 (Pinus tabulaefor- 
mis) 作 为 树 轮 样本 来 源 , 探 讨 油 松 生 长 对 气候 变化 
的 响应 ,并 基于 VS 模型 探讨 其 机 制 ,以 期 为 理解 气 
候 变 化 背景 下 该 区 域 油 松 对 气候 的 响应 机 制 和 评 
佑 气候 重建 的 可 靠 性 提供 依据 。 


1 研究 区 概况 和 方法 


1.1 研究 区 概况 

甘肃 省 天 水 麦 积 山 风景 名 胜 区 (34°07’~34°28'N， 
105°56'~106°10'E) 位 于 西 秦 岭 北 支 东 端 (图 1), 海 
$e 1450~2100 m, 紫 邻 黄土 高 原 , 北 邻 渭水 , 东 接 嘉 
陵 59 ,是 黄河 水 系 和 长 江水 系 的 分 水 岭 ; 气 候 上 属 
于 暖 温 带 湿 润 半 湿 润 气候 区 ,年 平均 气温 9~13 C, 
年 降水 量 420~660 mm, 多 集中 在 6 一 9 H ,无 霜 期 
170~220 d ,土壤 以 山地 标 培 为 主 ” ,植被 主要 属 暖 
温带 针 间 叶 混交 林 "… 。 
1.2 样 芯 采 集 及 年 表 建 立 

分 别 于 2021 年 3 月 和 5 月 两 次 在 麦 积 山北 坡 间 
山 瀑 布 (34.34°N,106.02°FE, 海 拔 1720 m) 设 立 采样 
点 ,在 胸 高 1.3 m 人 处 使 用 生长 锥 在 不 同方 位 采集 2 根 
FERS , 共 采 集 24 棵 油 松 ,获得 48 棵 样 芯 。 

在 实验 室 对 采集 的 树 芯 按 照 通用 的 树 轮 处 理 
方法 进行 预 处 理 , 利 用 LinTAB 年 轮 宽 度 测量 仪 在 
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0.01 mm 水 平 上 测定 轮 宽 , 并 使 用 COFFECHA 对 测 
量 结果 进行 检验 ,剔除 误差 较 大 的 样 芯 ,最 后 保留 
30 根 样 世 用 于 构建 年 表 。 利 用 ARSTAN 程序 的 负 
指数 函数 对 树 轮 宽度 数据 进行 去 趋势 和 标准 化 ,最 
终 得 到 差 值 年 表 (RES) 、 自 回归 年 表 (ARS) 及 标准 
化 年 表 (STD)。 标 准 化 年 表 主 要 排除 了 树木 自身 生 
长 趋势 造成 的 干扰 , 仅 保 留 大 部 分 高 .低频 环境 信 
号 ,选用 标准 化 年 表 进 行 研 究 ( 图 2)。 
1.3 气象 资料 

气象 数据 选用 的 是 离 采 样 点 最 近 的 麦 积 气象 
站 (34.57°N , 105.87°E , YFIR 1085.2 m) 的 逐日 观测 数 
据 , 数 据 来 源 于 资源 环境 科学 数据 中 心 (https://www. 
resdc.cn/Defaultaspx) 提 供 的 中 国 气 象 要 素 站 点 观 
测 逐 日 数据 集 。 根 据 麦 积 气 象 站 1980 一 2019 年 观 
测 资 料 ( 图 3) , 麦 积 站 多 年 平均 年 降水 量 523.13 
mm,8 月 最 多 (92.94 mm) ,12 月 最 少 (3.04 mm); 多 
年 平均 气温 11.45 % ,7 月 最 高 (23.33 C),1 月 最低 
(-1.71 ©), 

已 有 研究 指出 20 世 纪 90 年 代 西 秦岭 区 域 普遍 
发 生 了 气温 突变 中 通过 Mann-Kendall 检 验方 
法 对 麦 积 山 1980 一 2019 年 的 年 平均 气温 进行 分 析 
(图 4) ,1997 年 麦 积 山 气 温 发 生 了 突变 ,这 与 已 有 研 
究 结果 基本 一 致 ”“。Babst 等 "指出 树木 生长 与 气 
候 因 子 的 响应 模式 与 气温 有 关 , 因 此 ,本 文 将 1997 
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图 1 麦 积 山 采 样 点 及 气象 站 点 分 布 


Fig. 1 Distribution of tree-ring sampling site in Maijishan Mountain and meteorological station 
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图 2 树 轮 指数 与 样本 量 
Fig. 2 Tree-ring index and sample size 
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图 3 1980 一 2019 年 多 年 平均 月 均 温 和 月 降水 量 


Fig. 3 Monthly averaged temperature and precipitation 
during the period 1980 to 2019 


年 作为 研究 油 松 生 长 与 气候 因子 关系 的 时 间 分 界 
点 (1980 一 1997 年 ,1998 一 2019 年 )。 
1.4 VS 模型 
VS 模型 自 诞 生 以 来 已 经 开发 和 改进 多 个 版 本 ， 
本 文采 用 的 是 VS-GENN (http:/vs-genn.ru/) 提 供 最 
新 的 VS-oscilloscope, 其 与 过 去 版 本 在 基本 原理 上 一 
致 ，, 即 将 程序 分 为 4 个 模块 ,分 别 为 数据 输入 模 
块 、 基 本 模块 .用 来 估计 各 季节 细胞 数量 和 大 小 的 
模拟 形成 层 模块 和 数据 输出 模块 。VS-oscilloscope 
计算 树木 生长 速率 采用 的 公式 如 下 : 
Gr(t) = GrE() x min{GrT(,GrW(0} (1) 
式 中 : 6cr() 为 树木 年 轮 总 体 生 长 速率 ;CrE( 、Cr7TO)、 
Cr 人 分 别 为 光照 五 .气温 了 7 和 土壤 水 分 到 引发 的 逐 
日 生长 速率 。 
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图 4 1980 一 2019 年 麦 积 山 年 平均 气温 的 M-K 检 验 
Fig. 4 The M-K test of the annual average temperature in 
Maijishan Mountain from 1980 to 2019 


VS-oscilloscope 与 过 去 版 本 的 VS 模型 主要 差别 
是 改进 了 GrE 的 计算 公式 ,在 不 考虑 融 土 的 条 件 
下 ,采用 的 水 平面 上 每 日 接收 的 光照 的 计算 公式 
为 : 
E=I,, X(cosL X cosô X sinw + w X sin L X sin 8) Xxrx24/m 
(2) 
式 中 :天 是 太阳 常数 ;r 是 地 球 大 气 层 外 正常 人 射 的 
太阳 辐射 强度 与 太阳 常数 的 比值 ,无 量 纲 记 是 纬度 
(?);6 是 太阳 和 斜 角 (");o 是 日 落 时 角 ( 弧 度 )。 


2 结果 与 分 析 


2.1 年 表 统计 特征 分 析 

从 油 松 标准 化 年 表 统 计 参 数 表 中 可 以 看 出 ( 表 
1) ,年 表 有 较 高 的 平均 敏感 度 ,表明 油 松 对 气候 要 
素 变 化 敏感 ;一 阶 自 相 关系 数 较 低 ,表明 油 松 受 前 
一 年 气候 要 素 “ 滞 后 效应 ”的 影响 较 小 ; 样 芯 平均 、 
树 间 和 树 内 相关 系数 均 显 著 (P<0.05) , 信 噪 比较 高 ， 
第 一 主 成 分 为 0.377, 总体 代 表 性 大 于 90% ,这 些 都 
表明 年 表 具 有 较 高 的 可 信和 度 。 同 时 ,以 0.85 作为 总 
体 代 表 性 的 最 低 冰 值 …( 网 2) ,以 此 为 依据 建立 可 
靠 的 年 表 (1980 一 2019 年 )。 
2.2 树 轮 指数 与 气候 要 素 的 波动 特征 

将 1997 年 作为 研究 油 松 生 长 与 气候 响应 时 间 
的 分 界 点 (1980 一 1997 年 ,1998 一 2019 年 ) ,并 通过 
独立 样本 ;检验 比较 两 个 时 段 气候 因素 和 树 轮 指数 

异 ”。 由 麦 积 山 平均 气温 .降水 和 树 轮 指数 在 两 

个 时 段 的 箱 线 图 可 以 发 现 (图 5),1998 一 2019 年 的 
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表 1 年 表 统 计 特征 及 公共 区 间 分 析 
Tab. 1 Chronology statistics and the results of common interval analysis 
年 表 统计 项 数值 公共 区 间 分 析 数值 
年 表 人 研制 样本 量 30 公共 区 间 1993 一 2020 年 
序列 长 度 1947 一 2020 年 所 有 样 世 平 均 相 关系 数 0.321 
指数 均值 0.956 树 间 平均 相关 系数 0.315 
标准 差 0.209 树 内 平均 相关 系数 0.550 
平均 敏感 度 0.235 信 噪 比 14.173 
一 阶 自 相关 系数 0.234 总 体 代 表 性 0.934 
第 一 主 成 分 解释 量 0.377 
13.0 (a) 气温 900 - (b) 降水 1.8 『 (c) 树 轮 指数 
12.5 T 800 L is 。 极端 异常 值 
© 120 E 700 | 
z £ F 1.2 
a 11.5 = 600 L 8 7 
Br 110 j = 500 | 09 
10.5 400 上 i 
10.0 300 + 
1980—1997 1997—2019 1980—1997 1997—2019 1980—1997 1997—2019 
年 份 
图 5 1980 一 2019 年 树 轮 指数 与 气候 要 素 的 箱 线 图 
Fig. 5 Boxplots of tree-ring chronology and climatic variables in 1980-2019 
气温 显著 的 高 于 1980 一 1997 年 (P<0.05) ,降水 和 树 显著 性 水 平 


轮 指数 的 变化 则 不 显著 。 从 气候 要 素 上 ,1980 一 
1997 年 与 1998 一 2019 年 最 明显 的 差异 是 气温 的 显 
车 增加 。 

从 气温 与 降水 的 变化 趋势 上 可 以 看 出 (图 6)， 
除 1 月 10 月 11 月 外 ,其 他 各 月 平均 气温 和 全 年 平 
均 气 温 都 至 少 在 90% 的 显著 性 水 平 上 显著 增加 ; 而 
降水 仅 有 2 月 在 90% 的 显著 性 水 平 上 增加 。1980 一 
2019 年 麦 积 山 暧 十 化 趋势 显著 ,这 意味 着 树木 生长 
可 能 受到 更 强 的 水 分 胁迫 。 
2.3 树 轮 指数 对 气候 要 素 的 响应 

由 于 标准 化 年 表 的 一 阶 自 回归 系数 较 低 , 麦 积 
山 油 松 生长 受 前 一 年 气候 要 素 “ 浪 后 效应 ”影响 较 
小 。 因 此 ,选取 当年 逐 月 气温 、 降 水 与 标准 化 年 表 
进行 响应 分 析 。 将 1980 一 2019 年 .1980 一 1997 年 和 
1998 一 2019 年 3 个 时 间 段 的 气候 要 素 与 标准 化 年 表 
分 别 进行 相关 分 析 ( 图 7) ,1980 一 2019 年 油 松 生长 
与 5 月 气温 (P<0.01) .7 月 气温 (P<0.05) 呈 显著 的 负 
相关 ,与 6 月 降水 (P<0.1) .7 月 降水 (P<0.05) 呈 显著 
的 正 相 关 ,表明 5$ 一 7 月 高 温 导致 的 干旱 可 能 限制 了 


降水 


气温 


99% 增加 
F 95% 增加 
90% 增加 
不 显著 
90% 减少 
95% 减少 
123 45 6 7 8 910 11 12 WY 
月 份 99% 减少 
注 :WY 指 全 年 趋势 。 
图 6 1980 一 2019 年 气候 要 素 变化 趋势 
Fig.6 The trend of climatic variables change in 1980-2019 


油 松 的 生长 ,此 时 降水 使 土壤 湿润 程度 增加 可 以 组 
解 干旱 ,促进 油 松 生长 。1980 一 1997 年 油 松 生长 与 
5 月 气温 (P<0.01) .6 月 气温 (P<0.1) ,8 月 气温 (P< 
0.1) 旺 负 相 关 ,6 月 降水 (P<0.05) 呈 显著 正 相 关 , 整 
体 上 与 1980 一 2019 年 区 别 不 大 , 均 是 末 春 与 夏季 气 
温 导 致 的 干旱 对 油 松 生长 产生 了 限制 ,而 降水 有 助 
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图 7 油 松 标准 化 年 表 与 各 月 气候 要 素 的 相关 性 


Fig. 7 Correlations coefficients between he standardized chronology and monthly climatic variables 


于 缓解 干旱 ,促进 油 松 生 长 。1998 一 2019 年 的 响应 
模式 较 1980—1997 年 则 有 了 较 大 变化 ,5 月 .6 月 、8 
月 的 气温 和 6 月 降水 与 油 松 生长 的 关系 均 不 显著 ,7 
月 气温 、 降 水 与 油 松 生 长 的 关系 均 由 不 显著 变 为 显 
著 (P<0.05); 秋 冬季 节 的 气象 要 素 与 油 松 生长 表现 
出 了 相关 性 :10 月 气温 (P<0.1) 表 现 为 正 相关 ,12 月 
气温 (P<0.05) 、10 月 降水 (P<0.01) 表 现 为 负 相 关 。 
油 松 生 长 与 气候 的 响应 模式 在 1997 年 后 发 生 了 明 
显 改 变 。 

为 了 进一步 分 析 油 松 生长 与 气候 要 素 相关 性 
的 变化 ,选择 1997 年 前 后 表现 出 明显 变化 的 5 月 、7 


月 .10 月 、12 月 气温 和 6 月 、7 月 、10 月 降水 进行 10 a 
窗口 的 向 前 滑动 相关 分 析 ( 图 8) ,5 月 气温 与 油 松 生 
长 的 相关 性 不 断 下 降 ,10 月 12 月 气温 ,7 月 、10 月 
降水 的 相关 性 不 断 增加 ,而 6 月 降水 .7 月 气温 的 相 
关 性 则 表现 出 趋势 相似 的 震荡 。 这 表明 麦 积 山 油 
松 生 长 与 气候 响应 模式 变化 有 相关 性 下 降 、 相 关 性 
上 升 和 相关 性 震荡 3 种 模式 。 
2.4 VS 模型 模拟 

VS 模型 需要 输入 麦 积 山 逐 日 气温 降水 数据 . 采 
样 点 的 纬度 和 生长 参数 ,通过 调整 合适 的 生长 参数 
使 模型 模拟 结果 接近 实测 数据 路 。 生 长 参数 可 
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图 8 年 表 序 列 与 气候 要 素 的 向 前 滑动 相关 分 析 


Fig. 8 Correlation analysis of forward evolutionary intervals between chronology and climatic variables 
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通过 实测 资料 和 模型 模拟 调整 确定 ,模型 模拟 是 
根据 已 有 的 参数 范围 给 出 估计 的 参数 初始 值 , 根 
据 模 型 模拟 的 结果 与 实测 值 进行 对 比 ,不 断 调整 模 
型 参数 , 当 模拟 值 与 实测 值 一 致 性 够 高 时 认为 此 时 
的 各 个 参数 的 数值 基本 合理 。 由 于 一 些 生 长 参数 
对 树木 生长 的 作用 相互 持 抗 ,这 种 方法 获得 的 参数 
存在 错误 的 可 能 性 。 因 此 ,本 文 查阅 了 研究 区 的 
植被 ”土壤 “气候 ”等 资料 ,结合 实测 结果 和 
模型 模拟 共同 确定 合适 的 生长 参数 ( 表 2)。 
采用 表 2 中 列 出 的 参数 对 油 松 生长 进行 模拟 
(E19) ,1980 一 2019 年 模型 模拟 与 实测 值 的 相关 性 
为 0.59(P<0.05) ,1980 一 1997 年 和 1998 一 2019 年 也 
均 达 到 了 95% 显 车 性 水 平 ,说 明 模 拟 拟 合 程度 较 
好 ,可 以 真实 地 反映 麦 积 山 油 松 生长 情况 "。 
2.5 树 轮 指数 对 气候 要 素 的 响应 及 变化 的 解释 
FH VS 模型 模拟 的 逐日 平均 总 体 生 长 速率 (Cr) 
MERE (GrW) .气温 (Cr7) 及 光照 (GrE) 引 发 的 
生长 速率 可 以 看 出 (图 10) ,由 气温 引发 的 生长 速率 
E UEA ,在 150~250 日 (5 一 8 月 ) ,气温 对 油 松 生 
长 速率 的 贡献 下 降 , 这 是 由 于 气温 已 经 超过 VS 模型 
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Ko 8 月 气温 在 1980 一 1997 年 与 油 松 生 长 的 负 相 关 
原因 可 能 是 高 温 时 气温 超过 油 松 最 适 生 长 气温 上 
限 ,高 温 对 油 松 生 长 起 到 抑制 作用 ,推迟 了 和 气温 引 
发 的 生长 速率 的 第 二 个 峰值 所 致 。 

气温 引发 的 生长 速率 谷底 在 6 一 7 月 200 日 前 
后 。 由 于 干旱 缺 水 ,土壤 湿度 引发 的 生长 速率 也 在 
100 日 左右 开始 下 降 , 并 在 250 日 左右 下 降 到 最 低 
值 。 此 时 高 温 可 以 进一步 降低 气温 引发 的 生长 速 
率 , 加 剧 土壤 的 干旱 ,而 降水 可 以 缓解 干旱 ,提升 土 
壤 湿 度 引 发 的 生长 速率 。 这 解释 了 6 月 .7 月 油 松 生 
长 与 气温 的 负 相关 .与 降水 的 正 相关 。 

为 了 分 析 油 松 生长 过 程 的 变化 ,分 析 了 麦 积 山 
油 松 1980 一 1997 年 与 1998 一 2019 年 逐日 生长 速率 
的 差 值 和 1980 一 2019 年 生长 季 的 变化 (图 11), 
1998 一 2019 年 和 1980 一 1997 年 相 比 ,全 年 土壤 湿度 
引发 的 生长 速率 都 有 所 减少 ,反映 油 松 受 到 了 更 强 
的 水 分 胁迫 。 气 温 引 发 的 生长 速率 在 30~150 日 和 
270~330 H , , 即 2 一 $ 月 的 末 冬 和 春季 .9 一 11 月 的 秋 
季 有 所 增加 ;在 150~270 日, 即 5 一 9 月 的 春 末 和 夏 
季 有 所 下 降 ; 而 在 年 初 到 30 日 .330 日 到 年 末 变 化 不 


所 给 出 的 油 松 最 适 生长 气温 上 限 (23%)。 当 5 月 气 
温 过 高 时 ,就 有 可 能 使 气温 引发 的 生长 速率 的 峰值 
提前 ,对 油 松 生长 起 到 负面 作用 ,这 解释 了 为 何 
1980—2019 年 整体 上 5 月 气温 与 油 松 生长 呈 负 相 


大 。 这 些 变化 是 研究 区 变 暖 所 导致 ,2 一 5 月 .9 一 11 
月 没有 充足 热量 供 油 松 生长 ,因此 ,2 一 5 月 和 9 月 的 
增 温 使 油 松 生长 的 热量 供应 相对 更 加 充足 ,最 终 导 
致 生长 季 延 长 (P<0.05);5 一 9 月 的 气温 往往 已 经 超 


表 2 麦 积 山 油 松 生理 过 程 模拟 参数 


Tab. 2 Physiological process simulation parameters of Pinus tabulaeformis in Maiji Mountain 


参数 描述 值 
T AR BAR A/C 5 
To 最 适 生 长 气温 下 限 /%C 15 
Ts 最 适 生 长 气温 上 限 /%C 23 
Tı Ekim aC 33 
W, 初始 土壤 湿度 0.225 
Pix 饱和 土壤 的 最 大 日 降水 量 /mm 82 
1.( Droot) 根深 /mm 1600 
Cs(K,) 蒸腾 作用 第 一 系数 0.1175 
C3(K3) 蒸腾 作用 第 二 系数 0.185 
Ci(K) 降水 渗透 到 土壤 的 系数 (未 被 截流 的 系数 ) 0.74 
wW 生长 最 低 土 壤 含 水 量 (与 饱和 土壤 含水 量 比 ) 0 
LA 最 适 生长 土壤 含水 量 下 限 0.475 
Ws 最 适 生长 土壤 含水 量 上 限 0.725 
W, 植物 生长 最 高 湿度 (与 也 和 含水 量 比 ) 0.95 
CA(Ki) 土壤 排水 系数 0.001 
TT es) he OF RE KAER Y/C 10 
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Te -模拟 人 的 生长 速率 基本 一 致 , 总 体 生长 速率 差 值 与 气温 引 
发 的 生长 速率 差 值 的 相关 性 (0.66) 也 要 大 于 与 土壤 
ol 44 P 湿度 引发 的 生长 速率 差 值 的 相关 性 (0.45), 这 反映 

5 了 气温 引发 的 生长 速率 的 变化 是 导致 总 体 生长 速 

D 率 变化 的 主要 因素 。 但 总 体 上 生长 速率 总 体 变化 
os} 与 气温 引发 的 生长 速率 的 变化 也 有 区 别 ,主要 体现 
“él 在 冬 末 和 春季 的 增加 较 气温 引发 的 生长 速率 结 

得 更 早 (90 日 左右 ) ,秋季 的 增加 则 开始 较 晚 (290 日 
0 左右 ) ,这 反映 了 总 体 上 生长 速率 受到 土壤 湿度 引 


图 9 模拟 轮 宽 序列 和 实测 轮 宽 序列 对 比 
Fig. 9 Comparison of simulated and measured radial growth 


series 


过 VS 模型 所 给 出 的 油 松 最 适 生 长 气温 上 限 
(23 °C) ,此 时 继续 增 温 对 油 松 的 限制 更 强 ; 非 生长 
季 的 12 月 和 1 月 增 温 并 不 显著 ,因此 ,对 油 松 生长 
影响 不 大 。 总 体 生长 速率 的 变化 趋势 与 气温 引发 


发 的 生长 速率 减少 的 影响 , 增 温 导 致 的 干旱 抵消 了 
部 分 充分 热量 对 油 松 生长 的 促进 作用 。 

气候 变化 导致 的 油 松 全 年 生长 速率 的 变化 使 
油 松 生长 与 气候 响应 模式 发 生 了 改变 。5 月 气温 与 
树 轮 指数 负 相关 性 下 降 的 原因 可 能 是 2 一 5 月 变 暖 
导致 油 松 生长 的 热量 供应 更 充分 ,气温 引发 的 生长 
速率 更 高 ,抵消 了 5 月 高 温 的 负面 影响 。7 月 降水 与 
树 轮 指数 正 相关 性 的 增加 可 能 是 7 月 暖 干 化 使 高 温 


总 体 生长 速率 (Gr) 一 一 一 土壤 湿度 引发 的 生长 速率 (GrW) 
一 一 一 气温 引发 的 生长 速率 (Gr7) 光照 引发 的 生长 速率 (GrE) 
1.2 「 (a)1980—20194F 1.2 E (b)1980—19974F 1.2 F (c)1980—20194F 
1.0 1.0 
0.8 0.8 
K 
i 0.6 0.6 
MA 
H 0.4 0.4 
0.2 0.2 
上 n Aa oha NA o S SSS 
0 50 100 150 200 250 300 350 0 50 100 150 200 250 300 350 0 50 100 150 200 250 300 350 
日 数 日 数 日 数 


图 10 模拟 的 麦 积 山 油 松 生长 速率 


Fig. 10 Simulated growth rate of Pinus tabulaeformis in Maijishan Mountain 


(a) 生长 速率 差 值 


— Gr 


生长 速率 差 值 
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b) 生长 季 变化 


1980 1990 2000 2010 2020 


注 :SOS 为 生长 季 开 始 日 ;EOS 为 生长 季 结 束 日 ;DOS 为 生长 季 持 续 时 间 。 
图 11 麦 积 山 油 松 生 长 速率 差 值 及 生长 季 变 化 


Fig. 11 Growth rate difference and change of growing season of Pinus tabulaeformis in Maijishan Mountain 
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和 土壤 湿度 低 对 油 松 生长 的 胁迫 更 严重 ,气温 和 土 
二 湿度 引发 的 生长 速率 都 更 低 , 降 水 可 以 降低 气 
温 , 增 加 土壤 湿度 ,增加 气温 和 土壤 湿度 引发 的 生 
长 速率 ,最 终 促进 了 油 松 的 生长 。10 月 气温 与 树 轮 
指数 的 正 相 关 、 降 水 与 树 轮 指数 负 相 关 性 增加 的 原 
因 可 能 是 10 月 气温 引发 的 生长 速率 相对 增加 对 油 
松 生 长 起 到 促进 作用 ,促进 作用 超过 了 土壤 湿度 引 
发 的 生长 速率 相对 减少 对 油 松 生长 的 抑制 作用 , 导 
致 总 体 生长 速率 差 值 由 减少 转变 为 增加 ;此 时 更 高 
的 气温 可 以 增加 气温 引发 的 生长 速率 ,促进 油 松 生 
长 ;降水 伴随 的 降温 使 气温 引发 的 生长 速率 下 降 ， 
延缓 了 总 体 生长 速率 差 值 由 减 到 增 的 转换 ,对 油 松 
生长 产生 负面 影响 。12 月 植物 处 于 休眠 状态 ,气温 
很 难 对 油 松 生长 过 程 产生 影响 ,因此 ,12 月 气温 与 
树 轮 指 数 负 相关 性 的 增加 可 能 仅 是 相关 关系 而 非 
因果 关系 。6 月 降水 .7 月 气温 与 树 轮 指 数 相关 性 的 
震荡 变化 则 可 能 反映 了 研究 区 气候 存在 年 际 波动 
和 周期 性 变化 。 


3 讨论 


理解 树 轮 -气候 响应 关系 不 稳定 的 “分 离 效 应 ” 
和 更 本 质 的 树木 生长 的 气候 驱动 因子 的 变化 有 助 
于 预测 气候 变化 对 全 球 生 态 系 统 物 质 循环 和 能 量 
流动 的 影响 中 。 本 文 研究 了 天 水 麦 积 山 区 域 1980 一 
2019 年 油 松 生 长 对 气候 响应 的 变化 ,发 现 1997 年 前 
后 油 松 生长 与 气候 响应 关系 发 生 了 明显 改变 。 

树 轮 -气候 响应 关系 可 以 与 基于 树 轮 的 气候 重 
建 结果 相对 照 ,通过 验证 历史 时 期 树 轮 -气候 响应 
关系 的 稳定 性 以 评估 气候 重建 的 科学 性 中 。Chen 
等 Song 等 是 分别 基于 天 水 石门 山 采 集 的 油 松 样 
本 分 析 了 其 与 1952 一 2011 年 .1953 一 2008 年 气象 要 
素 的 响应 关系 , 均 发 现 与 5 月 .6 月 .7 月 的 气温 存在 
相关 关系 ,并 分 别 重建 了 天 水 1666—2008 年 5—6 
月 最 高 气温 和 1630 一 2011 年 5 一 7 月 平均 气温 。 这 
表明 5 月 .6 月 .7 月 气温 是 天 水 油 松 生长 的 主要 限 
制 因 子 。 本 研究 发 现 ,5 月 气温 与 油 松 生 长 的 相关 
性 在 1991 一 2000 年 后 才 有 比较 明显 的 下 降 ,1998 一 
2007 年 后 相关 性 才 下 降 到 95% 显 车 性 水 平 以 下 ,此 
前 两 者 相关 性 比较 稳定 。 因 此 , 树 轮 气候 重建 中 树 
轮 - 气 候 响应 关系 保持 稳定 的 假设 基本 成 立 , Chen 
等 9 Song 等 中 所 进行 的 气候 重建 比较 可 靠 。 
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Wu 等 ”基于 VS 模型 模拟 了 天 水 在 内 的 甘肃 东 
部 3 个 样 点 的 油 松 生 长 ,本 研究 中 VS 模型 参数 与 
Wu 等 ”选择 的 参数 主要 区 别 在 油 松 根深 及 相关 的 
土壤 湿度 参数 上 ,Wu 等 "选择 的 油 松 根系 深度 为 
500 mm, 本 研究 选择 的 油 松 根系 更 深 (1600 mm)。 
天 水 ”与 麦 积 山 生境 类 似 的 甘肃 武 都 五 凤 山 * 油 
松 生 理 的 研究 中 均 已 证 明 油 松 根深 可 以 达到 1500~ 
1600 mm ,更 深 的 根系 可 以 使 油 松 吸 收 深层 土壤 的 
水 分 ,对 油 松 的 生理 过 程 产生 影响 。 因 此 ,将 模型 
中 油 松 根系 深度 设置 为 1600 mm 可 以 更 好 的 反映 油 
松 的 实际 生长 情况 。 本 研究 也 对 一 些 土壤 湿度 参 
数 进 行 了 调整 。 史 江 峰 等 汪 刀 基于 VS 模型 对 贺 兰 
山 油 松 的 研究 中 ,也 对 油 松 根系 深度 和 土壤 湿度 参 
数 进行 过 类 似 的 修正 。 

基于 VS 模型 ,本 研究 发 现 气候 变化 使 油 松 全 年 
生长 速率 发 生 了 变化 ,进而 改变 了 油 松 生长 与 气候 
响应 的 关系 。 麦 积 山 变 暖 使 春秋 季 油 松 生 长 的 热 
量 供 应 充分 ,显著 地 提前 和 延 后 了 生长 季 开 始 日 和 
结束 日 ;同时 夏季 油 松 生长 受到 更 强 的 高 温和 干旱 
胁迫 。 生 长 季 开 始 和 结束 时 间 附 近 几 日 的 气温 、 降 
水 往往 对 油 松 生长 产生 较为 显著 的 影响 ,生长 季 变 
化 改变 了 各 月 气候 因子 与 油 松 生长 的 相关 性 ,突出 
表现 为 5 月 气温 与 树 轮 指数 负 相 关 性 的 下 降 、10 月 
气温 与 树 轮 指数 正 相 关 性 和 降水 与 树 轮 指数 负 相 
关 性 的 上 升 。 由 此 可 见 , 麦 积 山 油 松 生长 与 气候 响 
应 模式 的 变化 主要 是 由 当地 气候 变化 (特别 是 变 
暧 ) 导 致 ;如 果 变 上 暧 持 续 , 响 应 模式 可 能 发 生 进一步 
的 变化 。 考 虑 到 气候 变 暧 和 油 松 分 布 在 我 国 的 普 
遍 性 ,特别 是 一 些 区 域 油 松 生长 模式 与 研究 区 较为 
相近 ,类 似 的 油 松 生 长 与 气候 响应 模式 变化 在 其 
他 区 域 可 能 也 同样 存在 ”1。 

尽管 利用 VS 模型 探讨 树木 生长 过 程 的 研究 已 经 
非常 常见 ,但 是 模型 模拟 并 不 能 完全 代替 实际 观测 , 
也 不 能 保证 模型 模拟 的 所 有 生理 参数 都 合理 1。 
因此 ,如 果 要 进一步 研究 麦 积 山 油 松 生长 及 其 与 气 
候 因 子 的 响应 模式 ,需要 对 油 松 生长 开展 长 期 定位 
试验 和 连续 监测 ,利用 高 精度 实际 观测 数据 校正 模 
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背景 下 ,天 水 麦 积 山 区 域 油 松 树 轮 在 两 个 时 间 段 
(1980 一 1997 年 ,1998 一 2019 年 ) 对 气候 变化 的 响应 
差异 ,并 基于 VS 模型 对 其 变化 机 制 进行 了 探讨 , 主 
要 得 出 如 下 结论 : 

(1) 在 天 水 地 区 气候 显著 暖 干 化 的 背景 下 , 研 
究 区 油 松 与 气候 的 响应 关系 发 生 了 明显 改变 ,变化 
存在 相关 性 下 降 、 上 升 和 震荡 3 种 模式 ,说 明 研究 区 
油 松 生长 与 气候 的 关系 不 完全 稳定 。 麦 积 山 5 月 气 
温 与 油 松 生 长 的 相关 性 在 1991 一 2000 年 后 才 有 比 
较 明显 的 下 降 ,前 人 基于 天 水 石门 山 油 松树 轮 样 本 
重建 的 1666 一 2008 年 5 一 6 月 最 高 气温 和 1630 一 
2011 年 5 一 7 月 平均 气温 比较 可 靠 。 

(2) VS 模型 模拟 的 结果 表明 , 麦 积 山 油 松 生长 
与 气候 响应 模式 的 变化 主要 是 由 气候 变 暧 和 生长 
季 的 变化 引起 的 ,如 果 变 暖 持续 ,响应 模式 可 能 发 
生 进 一 步 的 变化 。 类 似 的 变化 在 有 油 松 生长 旦 生境 
类 似 ` 气 候 有 变 暖 趋势 的 其 他 区 域 可 能 同样 存在 。 

(3) 基于 VS 模型 模拟 的 麦 积 山 油 松 的 生长 过 
程 对 真实 生长 过 程 拟 合 效果 较 好 ,模拟 的 油 松 生长 
的 逐日 特征 和 各 因素 影响 也 较 好 地 解释 了 油 松 生 
长 与 气候 因子 间 响 应 关系 的 变化 机 制 。 因 此 ,可 以 
尝试 将 VS 模型 用 于 解释 其 他 区 域 的 树木 生长 与 气 
候 因子 间 响 应 关系 及 其 变化 。 
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The response and mechanism of Pinus tabulaeformis tree-ring width to climate 


change in Maijishan Mountain, Tianshui, China 


YAO Daijun, LIU Kang, HUI Yuxiang, WANG Kaixin 


(College of Urban and Environmental Science, Northwest University, Xi’ an 710127, Shaanxi, China) 


Abstract: In recent decades, an unusual reduction in forest radial growth and temperature sensitivity has been 
widely observed in the northern high latitudes. Several studies have also suggested that unstable and nonlinear 
relationships exist between tree growth and its climatic drivers at mid-latitudes and even globally. However, this 
relationship remains poorly understood, particularly in the mid-latitudes. The climate response and growth process 
before and after the temperature abrupt change were investigated using dendroclimatological methods and the 
Vaganov- Shashkin model in order to research the variations and mechanisms of response of radial growth of 
Pinus tabulaeformis to climate factors in Maijishan Mountain, Tianshui, during 1980- 2019. The findings 
revealed that: (1) According to the Mann-Kendall test, temperature changed abruptly in 1997 and then increased 
significantly. During 1980-1997 and 1998-2019, there were three patterns of response in radial growth to climate 
factors revealed in the variation in correlation with climatic factors: decline (Temperature in May), increase 
(Temperature and precipitation in October, precipitation in July and temperature in December), and fluctuation 
(Precipitation in June and temperature in July). (2) The Vaganov- Shashkin model estimated chronologies are 
significantly correlated with the measured tree-ring chronologies before and after the abrupt change (P<0.05). 
According to simulation results, significant climate warming resulted in more abundant heat supply for growth in 
spring and autumn, thereby significantly extending the growing season and potentially causing a change in 
response patterns in May and October. Furthermore, summer high temperatures and water scarcity impacted 
growth, potentially leading to a change in response patterns in July. (3) The change in response patterns to climate 
factors is primarily caused by climate warming and the lengthening of the growing season. If the climate 
continues to warm, the response patterns are expected to change, even more, and similar behavior may exist in 
other areas where the species is found. 


Keywords: tree-ring; climate change; Vaganov-Shashkin model; Pinus tabulaeformis; Maijishan Mountain 


